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Second  oraer  rigid  flapping  response  statistics  of  lifting 
rotor  blades  are  obtained  by  efficient  algorithms  developed  in 
this  paper.  These  statistics  enable  us  to  analyze  the  effect 
of  a finite  load  correlation  length  on  the  blade  response.  For 
our- particular  loading  (typical  of  turbulence  excitations)  and 
for  an  advance  ratio  small  compared  to  unity,  this  effect  can 
be  concisely  expressed  in  terms  of  an  amplitude  factor  and  a 
phase  factor . In  the  case  of  a blade  excited  by  a random 
vertical  inflow  with  a spanwise  correlation  length  of  the  order 
of  the  blade  length,  the  amplitude  factor  shows  that  there  may 
be  as  much  as  a 40%  error  in  a solution  which  assumes  the  irifloi 
is  spatially  uniform.  The  analytical  development  leading  to  th< 
algorithms  also  illustrates  how  a spatial  correlation  method 
for  general  linear  PDb  with  random  forcing  previously  formu- 
lated by  the  author  may  be  used  in  conjunction  with  a Ritz- 
Galerkin  procedure. 
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Abstract 


beconu  oruer  rigid  flapping  response  statistics  of  lifting 
rotor  blades  are  obtained  by  efficient  algorithms  developed 
in  this  paper.  These  statistics  enable  us  to  analyze  the 
effect  of  a finite  load  correlation  length  on  the  blade  response. 
For  our  particular  loading  (typical  of  turbulence  excitations) 
and  for  an  advance  ratio  small  compared  to  unity,  this  effect 
can  be  concisely  expressed  in  terms  of  an  amplitude  factor 
and  a phase  factojr.  In  the  case  of  a blade  excited  by  a 
ranaom  vertical  inflow  with  a spanwise  correlation  length  of 
the  order  of  the  blade  length,  the  amplitude  factor  shows  that 
there  may  be  as  much  as  a 40%  error  in  a solution  which  assumes 
the  inflow  is  spatially  uniform.  The  analytical  development 
leading  to  the  algorithms  also  illustrates  how  a spatial 
correlation  method  for  general  linear  PDE  with  random  forcing 
previously  formulated  by  the  author  may  be  used  in  conjunction 
with  a Ritz-Galerkin  procedure. 
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Nomenclature 


dimensionless  transverse  displacement  of  the  blade 
blaue  length 
Lock  number 
= y/6 

uniform  rotating  speed  of  the  blade 

constant  forward  velocity  of  vehicle-blade  system 

advance  ratio  (=  V^/USL) 

4 2 

dimensionless  stiffness  factor  (=  LI/ml  ft  ) 
linear  mass  density  of  the  blade 
inflow  ratio 

dimensionless  correlation  time  (see  eq.(3)) 

spatial  part  of  the  autocorrelation  of  the  loading 
(see  eq . ( 3) ) 

a temporally  delta  correlated  random  process 
(see  eq . (4) ) 

spatial  correlation  functions  (see  eq.(6)) 
autocorrelation  of  w(x,t) 

dimensionless  correlation  length  of  loading 
(see  eq.  (28) ) 

a positive  constant  (see  eq.(28)) 
see  eq.  (9)  , (10)  , (12)  and  (13) 

flapping  angle  of  the  blade 

see  eq.  (18)  (meansquare  flapping  angle,  etc.) 

the  spring  force  and  damping  parameter  (see  eq.(20)) 

see  eq.(19),  (21)  and  (22) 
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see  eq.(23)  and  (24) 

steady  state  solution  of  p and  q 

see  eq. (32) 

steady  state  solution  of  U , S and  V 

amplitude  and  phase  factor  (see  eq.(32)  and  (47)) 
in  the  steady  state  response  to  random  inflow 

see  eq . (31) 

0 and  V with  e = 0 

( ) for  a random  change  in  collective  pitch  0(x,t) 
= < A (x,t )w  (y , t)  > 
see  eq. (57) 

= <<j>  (T  ) <Mt  ' ) > 
see  eq . (60) 
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1.  Introduction 


The  dynamics  of  flexible  lifting  rotor  blades  in  forward 
flight  (Fiy.l)  is  complicated  by  the  fact  that  the  aerodynamic  lift 
acting  on  the  blade  changes  significantly  in  the  course  of  each 
blade  revolution.  Even  an  analysis  of  the  small  amplitude  motion 
of  such  a structure  must  cope  with  problems  such  as  parametric 
excitation  associated  with  the  periodically  time-varying  system 
parameters  which  characterize  the  aerodynamic  damping  and  spring 
force  effects.  In  one  model  for  the  forced  small  transverse 
vibration  of  a single  blade,  the  dimensionless  transverse  displace- 
ment w(x,t)  (normalized  by  the  blade  length  £)  is  governed  by 
the  rather  formidable  (dimensionless)  partial  differential 
equation  [1,2,9] 

Wtt  + Y0lX+PsinTK  + Lxt  ^ = f(x,x)  (0<x<l,  x>  0 ) (1) 

with 

Lxt[  ] = ^ ]xxxx  - Jxx  + (x+Y0MCOST|x+ysinx|) [ Jx 

(2) 

where  Y = 6Yq  is  the  Lock  number  characterizing  the  aerodynamic 
effect,  y is  the  advance  ratio  (the  ratio  of  the  forward  speed 
of  the  vehicle,  V^,  to  the  rotating  speed  at  the  blade  tip,  fi£) , 
x is  the  distance  from  the  axis  of  rotation  along  the  blade  span 
normalized  by  the  blade  length,  and  x/Sl  is  the  real  time.  The 

4 

effective  bending  stiffness  factor  of  the  blade,  £ , is  related 

to  the  bending  stiffness  of  the  uniform  blade,  El,  by  the  relation 
4 4 2 

C = EI/m£  (2  , where  m is  the  lij-iear  mass  density  of  the  blade. 


When  the  source  of  external  excitation  is  a vertical  inflow,  we 
have  f ( x , t ) = Yq  | | X (x,t  ) where  X is  the  so-called 

inflow  ratio.  The  temporally  periodic  coefficients  in  the  PUE  (1) 
give  rise  to  the  possibility  of  parametric  excitation  and  dynamic 
instability  (see  [1]  and  references  therein). 

Aside  from  the  various  stability  analyses,  there  is  also 
the  problem  of  the  effect  of  random  air  and  rotor  generated 
turbulence  on  the  structural  integrity  of  the  blade.  In  an 
effort  to  understand  this  aspect  of  the  rotor  blade  problem, 
several  recent  papers  studied  the  stochastic  blade  response  to 
a (zero  mean)  random  inflow  with  known  statistics  (see  [2]  and 
references  given  therein) . because  of  the  time-varying  coefficients 
in  (1),  the  steady  state  response  process  w(x,t)  will  be 
temporally  nonstationary  even  if  A(x,x)  is  stationary.  In  spite 
of  the  substantial  reduction  (by  at  least  an  order  of  magnitude) 
in  machine  computation  made  possible  by  a new  method  of  solution 
developed  in  [3]  and  used  in  [2],  i"  is  still  rather  expensive 
to  generate  useful  information  on  the  stochastic  properties  of 
the  nonstationary  steady  state  response  of  the  flexible  blade 
for  design  purpose. 

If  the  inflow  is  uniform  along  the  blade  span  so  that  X(x,x) 
is  independent  of  x,  one  may  expect  that  the  dominant  motion  of 
a blade  hinged  at  the  axis  of  rotation  is  in  the  form  of  rigid 
flapping.  A solution  of  the  stochastic  forced  transverse  vibra- 
tion problem  for  a spanwise  uniform  loading  based  on  a rigid  flap- 
ping blade  model  is  considerably  simpler  than  a more  general 
flexible  blade  analysis  as  far  as  the  amount  of  required  machine 
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computation  is  concerned  (see  [4]  and  [5]).  Inasmuch  as  the 
stochastic  loadings  experienced  by  rotor  blades  are  often  random 
functions  of  both  space  and  time,  a rigid  flapping  solution  for 
(zero  mean)  spanwise  correlated  random  loads  of  comparable 
simplicity  should  be  of  interest  (as  pointed  out  in  [4]).  Such 

4 

a solution  and  the  string  solution  (q  = 0)  of  [2]  together 
delimit  the  range  of  the  solution  for  any  flexible  blade  with 

4 

finite  bending  stiffness  (0  < £ < »)  . With  the  help  of  the 

spatial  correlation  method  of  [3],  we  can  now  formulate  an 
efficient  computational  procedure  to  obtain  such  a rigid  flapping 
solution  for  a spanwise  correlated  random  excitation  containing 
as  a special  case  the  solution  of  [4]  and  [5]  for  spanwise  uniform 
inflows.  While  our  results  constitute  a step  toward  a better 


understanding  of  rotor  blade  behavior  under  random  excitation, 
the  analytical  development  leading  to  these  results  also  illus- 
trates how  the  general  spatial  correlation  method  may  be  used  in 
conjunction  with  a Ritz-Galerkin  procedure. 

For  the  purpose  of  illustrating  our  method  of  solution,  we 
take  X(x,t)  to  be  of  zero  mean  and  exponentially  correlated 
in  time  with  an  autocorrelation  function 

““ I T 2_T 1 1 

<X  (x2,t2)A (x1,x1)>  = e Rg(x2,x1)  (3) 

where  <...>  is  the  ensemble  averaging  operation,  a is  a known 

positive  constant  and  Rs(x2,x^)  = Rgfx^.x^  is  a given  function. 

2 

Rs  (x1,x2)  was  taken  to  be  a positive  constant  a for  the  case 
of  a random  inflow  due  to  high  altitude  air  turbulence  in  [4]. 
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Since  equation  (1)  is  linear  (so  are  the  associated  initial  and 
boundary  conditions),  w(x,t)  is  also  of  zero  mean  and  we  can  there- 
fore concentrate  on  the  second  order  response  statistics  of  w(x,t) 
characterized  by  the  autocorrelation  function  R (x2  ,x  2 ’>  xy  > T 
= <w  (x2  , x2 ) w (x^ , x ^ ) > . To  determine  R (x2 ,x 2 ;x^ ,x ^) , we  con- 

sider A(x,x)  to  be  the  steady  state  stationary  response  to  a 
temporally  uncorrelated  random  excitation  n(x,x)  of  a dynamical 
system  characterized  by  the  first  order  ODD  (see  [2,3]) 

A^  + aA  = /2a  n(x,x)  (4) 

where  <n (x2 ,x 2 ) n (x^ ,x ^) > = R^ (x2 ,x^ ) 6 (x 2~x ^ ) . It  is  not  difficult 
to  verify  that  the  autocorrelation  function  of  the  steady  state 
solution  of  (4)  is  as  given  by  the  right  hand  side  of  (3) [2]. 
Furthermore,  it  can  be  shown  [2]  that 

<n(y,x  ' )w(x,x)>  = <n  (y  ,x  ' ) w^  (x  ,x ) > = 0 (5) 

for  all  x ' >_  x > 0 and  0 <_  x,y  <_  1.  The  numerical  results  to 
be  given  in  this  paper  will  be  for  the  special  case  R (x,y) 

= o exp  (~e  I x-y  | ) where  o > 0 and  c >_  0 are  given  constants. 

The  analytical  and  numerical  results  for  the  rigid  flapping 
solution  obtained  with  the  help  of  the  above  device  allow  us  to 
study  the  effect  of  a finite  load-correlation  length  characterized 
by  the  dimensionless  number  e (with  the  correlation  length  equal 
to  Jt/e)  on  the  second  order  statistics  of  the  (zero  mean)  flapping. 

3 

blade  response.  In  the  low  advance  ratio  range,  u <<  1,  a pertur- 
bation solution  shows  that  the  effect  of  the  correlation  length  may 
be  completely  described  by  an  amplitude  factor  p and  a phase  factor 
v;  both  factors  are  simple  functions  of  e . In  the  case  of  a random 
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vertical  inflow  with  a correlation  length  of  the  order  of  the  blade 
length  (e  = 0(1)),  we  see  from  the  expression  for  p(e)  that  the 
discrepancy  between  our  solution  and  one  ignoring  the  finite  spatial 
load-correlation  (as  in  [4]  and  [5])  may  be  as  much  as  40%  of  the 
former.  In  the  high  advance  ratio  range,  an  efficient  numerical 
solution  procedure  is  formulated  for  the  second  order  statistics 
of  the  periodic  steady  state  flapping  blade  response.  The  numeri- 
cal solution  obtained  by  this  efficient  procedure  shows  that  the 
effect  of  a finite  load-correlation  length  (0  < e < ®)  is 
qualitatively  similar  to  that  described  by  the  amplitude  and 
phase  factor  for  the  low  aavance  ratio  case. 

2 . Spatial  Correlation  Functions  for  Flexible  Blade  Response 

The  essential  feature  of  the  spatial  correlation  method  for  the 
second  order  response  statistics  proposed  in  [3]  and  used  in  [2]  and 
(6)  is  the  formulation  of  a nonstochastic  mixed  initial-boundary 
value  problem  for  the  four  unknown  spatial  correlation  functions  of 
the  response  process  w(x,x): 

u (x,y , t)  = <w  (x , x ) w (y , x ) > , s(x,y,x)  = <w  (x , t ) w (y  , t ) > 

(6) 

t (x ,y , t ) = <wT (x,x)w(y,x) >,  v(x,y,x)  = <wt (x,x)wt (y,x) 

for  all  0 <_  x,y  1 and  x ^ 0.  Note  that  these  spatial  correlation 
functions  contain  the  meansquare  response  properties  as  special  cases 
(when  y=x) . As  we  shall  see,  they  also  serve  as  the  initial  condi- 
tions for  a nonstochastic  mixed  initial  boundary  value  problem  for 
the  determination  of  the  autocorrelation  function  R (x2 , X£ ;x^ , x ^) , 

(section  7) . 

To  obtain  an  appropriate  set  of  equations  for  u,  s,  t and  v, 

^ J 
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we  observe  that 

UT  ~ <wx (x,x)w(y,x)>  + <w (x,t)wt (y,T) > = t + s (7) 

and 

tT  = v(x,y,T)  + <wtt (x,T)w(y,x)>  (8) 

where  we  have  made  use  of  the  fact  that,  within  the  framework  of 
meansquare  convergence,  differentiation  commutes  with  the  ensemble 
averaging  operation.  We  now  use  equation  (1)  to  eliminate  w 

TT 

from  (3)  so  that 

fcx  = v " Lxt[u]  “ Y0 |x+usinx | t + p(x,y,T)  (9) 

where 


f>  (x,y , x)  = Y0|x+psinT|<A  (x,T)w(y,T)> 

i . , , (10) 

= Y q 1 x+y  smx  | p (x , y , t ) 

Interchange  the  role  of  x and  y and  we  have  also 

sx  = v ” Ly T ( u ^ “ Y 0 | y+ysinx | s + p(y,x,x)  (11) 

Finally,  similar  manipulations  applied  to  the  expression  for 
vT  give 

VT  = “ LxtIs1  " Lyx  " Yq ( | x+y sinx | + | y+ysinx | ) v + §(x,y,x) 

(12) 


where 


^(x,y,x)  = Y0|x+usinx |q (x,y,x)  + yq |y+ysinx |q (y ,x,x)  (13) 

with  q (x ,y , x ) = <\ (x,x)wt (y,x)  >. 


-7- 


Equations  (7) , (9) , (11)  and  (12)  are  to  be  satisfied  in  the 

interior  of  the  semi-infinite  unit  square  column  (0<x,y<l),  x>0) 
in  the  x,y,x-space.  On  the  base  square  of  the  column,  x = 0,  we 
have  from  the  condition  of  no  initial  transverse  motion: 


u (x, y , 0 ) = s (x,y , 0)  = t(x,y,0)  = v(x,y,0)  = 0 


(14) 


We  will  not  be  concerned  with  the  appropriate  boundary  conditions 
on  the  four  walls  of  the  column  (given  in  [6])  as  they  do  not 
enter  into  our  analysis  of  blade  flapping. 

The  four  equations  (7) , (9) , (11)  and  (12)  contain  six  unknowns 

since  p(xfy,x)  and  q(x,y,x)  involve  the  unknown  w(x,x).  We 
need  two  more  equations  to  complete  the  system.  To  get  these,  we 
observe  that 

PT(x>y/T)  = <Ax (x,x)w(y,x) > + <A (x,x)wx (y,x)> 

(15) 

= - ap(x,y,x)  + q (x,y , x ) 

and 

qx(x,y,T)  = <Ax  (x,x)wT  (y,x)>  + <A  (x,  x ) w (y , x)  > 

= - (a+Ag | y+ysinx | ) q (x,y , x ) - L [p(x,y,x)J 

+ Y0 |y+psinx |RS (x,y)  (16) 


where  we  have  made  use  of  the  PDE  (1)  to  eliminate  w^^,  the  ODE  (4) 
to  eliminate  A^,  and  the  conditions  (5)  to  simplify  the  resulting 
equations.  The  initial  conditions 

p (x ,y , 0>  = q (x,y , 0)  = 0 (0<x,y<l) 


(17) 
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supplementing  (15)  and  (16)  follow  from  the  fact  that  the  blade 
experiences  no  transverse  (out  of  the  rotor  plane)  motion  up  to 
some  reference  time  x = 0.  Again,  we  need  not  give  here  the  specific 
boundary  conditions  for  p and  q at  y = 0 and  y = 1,  other  than 
noting  the  fact  that  they  involve  only  two  unknowns,  p and  q.  As 
such,  we  can  first  solve  (15)  and  (16)  for  p and  q in  the  y,x-space 
with  x as  a parameter,  and  then  use  the  result  in  (7),  (9),  (11) 

and  (12)  for  the  determination  of  the  other  four  unknowns.  We  note 
also  that  the  spatial  correlation  of  the  loading,  characterized 
by  R^(x,y),  enters  into  the  analysis  explicitly  only  through  its 
appearance  on  the  right  side  of  (16)  . 

3 . The  Rigid  Flapping  Motion 

We  now  introduce  the  rigid  flapping  assumption  by  taking 
w(x,t)  = x<j>(t),  so  that 

u(x,y,T)  = xyU(x),  s(x,y,x)  = xyS(x) 

(18) 

t (x,y , x ) = xyT(x),  v(x,y,x)  = xyV(x) 

where  U(x)  = <<{)2(x)>,  etc.,  and  equations  (7)  , (9),  (11)  and  (12) 
become  four  ODE: 

U = T + S 

T = V - [w2+k  (x ) ) U - c (X ) T + P(X) 

S = V - [w2+k (x ) ] U - c(X)S  + P(x)  (19) 

V = - [ui2+k  (X)  ] (S+T)  - 2c  (X ) V + 2Q(X) 


where 


k ( x ) = 3yny cost  jx+ysinx|xdx 


'.C- 


x+ysinx I x dx 


{P(t),Q(t)}  = 3y  x|x+ysinx | {p(x,x),q(x,T) }dx 


(p(x,T) ,q(x,x) } = 3 [ y{p (x,y ,x) ,q (x,y,x) }dy  (22) 

> 0 

2 

The  constant  u is  equal  to  1 if  the  blade  is  hinged  at  the 
blade  root  and  is  greater  than  unity  if  there  is  an  elastic  root 
restraint. 

The  quantities  p(x,x)  and  q(x,x)  are  determined  by 

= - ap  + q,  qT  = -[a+c(x)]q  - [w2+k(x)]p  + r(x,x)  (23) 
p (x,  0)  = q (x,  0)  = 0 (24  ) 


where 


r (x,x ) = 3y 


°l!yl 


y+ysinx |Rs(x,y)dy 


Equations  (23)  are  obtained  from  (15)  and  (16)  by  multiplying 
through  by  3yQy  and  integrating  over  the  interval  (0,1). 

The  general  procedure  is  to  solve  the  initial  value  problem 
(23)  and  (24)  with  x as  a parameter.  The  results  are  to  be 
used  in  the  integrals  on  the  right  side  of  equations  (21)  and  the 
integrals  evaluated  to  give  P(x)  and  Q(x).  Having  P and  Q, 
we  can  then  solve  the  four  equations  (19)  subject  to  the  initial 
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conditions 

U (0)  = S(0)  = T ( 0 ) = V (0)  = 0 (26) 

which  follow  from  (14)  . We  note  however  that  the  second  and 
third  equation  of  (19)  together  with  S(0)  = T ( 0 ) = 0 imply 
S(x)  = T(t)  for  all  t (consistent  with  the  fact  S = <$#>  = T)  . 
Therefore,  the  system  (19)  is  effectively  a system  of  three  equations 

U = 2S,  S = V - (oj2+k  (x ) ] U - c(t)S  + P 

, (27) 

V = -2[w  +k(x)]S  - 2c  (t ) V + 2Q 

The  dampiny  coefficients  c(t)  and  the  supplementary  spring 

rate  k(x)  due  to  the  aerodynamic  lift  have  been  calculated  in 

[9] . In  the  case  where  X is  independent  of  x,  we  have 
2 

Rg(x,y)  =o  (a  positive  constant);  the  corresponding 
r(x,x)  = r(x)  reduces  to  the  envelope  function  for  the  inflow 
ratio  term  given  in  [9] . 

4 . Exponential  Correlation  in  Space—  Hovering 

In  the  remaining  sections  of  this  paper,  we  restrict  ourselves 
to  a special  class  of  load  spatial  correlation  function 

Rs(x,y)  = o2e"elx"yl  (28) 

2 

where  a >0  and  e>_0  are  known  constants.  We  will  be  interested 
in  how  the  rigid  blade  solutions  obtained  in  [4],  [5]  and  [7]  for 
a spatially  pniform  random  excitation  are  modified  by  a finite 
load  correlation  length.  In  this  section,  we  consider  first  the 
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simpler  case  of  hovering. 

With  p = 0,  equations  (23)  become 


p = -op  + q-  qT  = -(a  + J)q  - o p + r 


r = 3Yno2  y2e-elX-yldy 

0 (30) 

= l^[2x2+e2x4-x2e"CX-x2eC  (x_1)  (l+e+  \-) ] 
e J 

Since  r is  independent  of  x,  the  steady  state  solution  of  (29), 
denoted  by  pg  and  qp , is  also  independent  of  x and  can  be 
obtained  simply  by  setting  pt  = qT  = 0.  The  resulting  algebraic 
equations  give 

p = r(x)/A  , q = ap  , A = (J2  + a2  + ay/8  (31) 

s s s 


Correspondingly,  we  have  from  (21) 


1 2 2 

p = 3yn  [ x2p  (x)d;:  = p(e) 

s '0  JQ  3GA  (32) 

7?  -f  f2  f 3 e4  e5 

Qs  = aPs  , p(e)  = -e  [e  (l+e+  —)  ~ (1_  T + "8~  " 20)  1 

which  are  also  independent  of  x.  It  follows  that  the  steady  state 
solution  of  (27),  denoted  by  U , S and  V , is  also  independent 
of  x.  By  setting  = 0 , we  have  immediately  from  (27) 

V - T PS  - p<£)  £ V(E) 

2  (33) 

U = Unp(e)  5 XSLJ8 p (e ) , i = 0 
0 36Aw 


k 
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Note  that  Vg  and  Uq  are  independent  of  e and  are  in  fact 
the  meansquare  velocity  and  displacement  known  for  the  case 
of  a spanwise  uniform  random  inflow  witn  the  same  exponential  time- 
correlation  [4,5,7].  The  factor  p(e)  in  (32)  and  (33)  may  there- 
fore be  thought  of  as  an  amplitude  factor  associated  with  a finite 
spanwise  correlation  length  of  the  inflow.  It  is  not  difficult  to 
verify  that  p(e)  -*•  1 and  p'  (e)  <0  as  e -*■  0 so  that  V/V^ 
and  U/Uq  decrease  with  increasing  e for  small  e.  A small  but 
positive  e means  a finite  correlation  length  which  is  long  compared 
to  the  blade  length.  On  the  other  hand,  p(e)  -*■  0 and  i ep(e)  9/5 
as  e -*■  °°,  so  that  the  solution  tends  to  that  of  a spatially  delta- 
correlated  inflow.  The  variation  of  the  amplitude  factor  P(£)  over 
the  whole  range  of  c is  given  in  Figure  2 where  we  have  plotted 
ep(e)/2  for  all  t > 2 in  order  to  compare  with  the  limiting  case 
of  spanwise  delta-correlated  inflow.  The  plot  shows  that  p is 
a monotone  decreasing  function  as  e increases.  Therefore,  the 
meansquare  flapping  displacement  and  velocity  decrease  with  decreas- 
ing spanwise  correlation  length  of  the  particular  class  of  inflows. 

For  blades  in  hover,  the  problem  with  a random  inflow  excita- 
tion is  a rather  artificial  one;  an  excitation  due  to  a randomly 
changing  collective  pitch  angle  6(x,x)  is  more  appropriate.  For 
this  case,  we  have  f(x,x)  = y ^ | x+psinx | 6 (x,x ) . if  8(x,x)  is 
exponentially  correlated  both  in  space  and  time  (as  given  by  (3)  and 
(28) ) , similar  calculations  for  p = 0 give 

^ °e(E>'  °e  = 127ff¥Li°0<c)'  §e  = 0 <34> 

64Aw 


where 
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,K1  4 5 6 7 

/^\  _ 1152  f/i  C _i_£  C j_  C \ 

P0(e)  -g-  1(1  “ 2T  To  72  252' 

e 

2 3 (35) 

-e"C(l+  c + — + —)] 

and  where  A is  as  given  in  (31) . The  variation  of  with  e is 

also  shown  in  Figure  2.  With  p q ( c ) ->  1 as  e -*■  0,  the  results  in 

(34)  tend  to  those  for  a spanwise  uniform  0 obtained  in  [7] . As 

e -*•  °°,  Pq(c)  tends  to  zero  while  -jep0(e)  tends  to  16/7  corresponding 

2 

to  the  case  R (x,y)  = a 6 (x-y)  . For  finite  values  of  e»  pQ(e) 
is  again  a monotone  decreasing  function  as  e increases.  Therefore, 
the  meansquare  flapping  displacement  and  velocity  are  also  reduced 
by  a shortening  of  the  spanwise  correlation  length  of  this  particular 
class  of  0 (x,t  ) . 

Before  leaving  the  hovering  case,  it  should  be  noted  that  the 
quantities  Vq  and  Uq  (for  both  kinds  of  random  excitations  consiuer- 

ed)  are  monotone  increasing  functions  of  y for  all  y 0 and  for  all 

2 2 
positive  values  of  uj  and  a.  In  the  realistic  range  of  y and  to  , 

Vq  increases  almost  linearly  with  y while  Uq  increases  quadratically 

with  y for  a = 0 (1)  and  is  nearly  linear  in  y only  for  broad  band 

excitations  (a  >>  1) . 

5 . Forward  Flight  at  Moderate  Advance  Ratios 
While  an  exact  elementary  solution  of  our  problem  was  obtained 
in  section  (4)  for  the  hover  case  (y=0) , the  same  is  not  possible 
for  the  forward  flight  case  (y>0)  . To  gain  some  insight  into  the 
effect  of  a spanwise  correlation  of  the  inflow,  we  restrict  ourselves 

3 

in  this  section  to  the  low  advance  ratio  range  so  that  y <<  1.  In 
this  range,  we  expect  that  the  contribution  of  the  reverse  flow  ef- 
fect can  be  neglected  (see  [4.8])  so  that 


c (t ) = 3y, 


3 2 Y Y • — 

: +x  ysinx)dx  “ g + ^ usin  H cn(t) 


(CONTINUED) 
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r1  2 

k(x)  = 3yQycosx  (x  +xysinx)dx 
> 0 


= ^-pcosi  + Xy2sin2x  = k (t) 
do  n 


r (x,  t ) = 3y 


2 f1 
°°  Jo 


(y  +yysinx)e 


-c |x-y 1 


dy 


= a y [rQ (x)+r1 (x)psinx]  5 rn(x,x) 


where 


(36) 


rQ(x)  = c 3[2+e2x2-e“cx~e  c(1_x)(l+e+  |e2)] 
rx(x)  = je~2 [2ex+e“cx-e"c (1-x)  (1+e) ) 


(37) 


The  form  of  cn(x),  k (x)  and  rn(x,x)  suggests  that 

a steady  state  solution  of  (15)  and  (16)  in  powers  of  y is 

★ 

possible  when  y < 1.  Upon  writing 


= I (Pn  (x,x)  ,qn  (x, x)  }yn  (38) 

n=0  n n 


the  coefficients  p^  and  q^  (which  are  independent  of  y)  are 
evidently  the  particular  solutions  of  the  following  sequence  of 
ODE: 


PO  + 2C0P0  + Apo 

Pl  + 2cqP;  + APl 
ql  = pi  + apx 


Yr0(x)  , qQ 
yr^ (x) sinx  - 


= Pq  + apQ  (39a) 

|-sinxpQ  - ^(asinx+cosx)  pQ 

(39b) 


A perturbation  solution  of  the  initial  value  problem  (15) -(17)  it- 
self can  be  obtained  without  difficulty.  But  we  are  not  interested 
here  in  the  transient  part  of  the  meansquare  response  properties. 


*,.*-**>. 
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where  dots  indicate  differentiation  with  respect  to  x,  A is 
as  defined  in  (31)  and  Cq  = a + y/lG.  It  is  a straightforward 
matter  to  obtain  these  particular  solutions  since  the  ODU  involved 
are  with  constant  coefficients. 

The  steady  state  perturbation  solutions  (38)  are  then  inserted 
into  (21)  with  the  reverse  flow  effect  neglected.  Upon  carrying 
out  the  integration,  we  get 


P(t)  = — p (s)  {1+y  [PsQ+Pslv  (e)  ] sinx 


+ P lpco+pciv  COST  + °(P  )) 


Q(x)  = P (e)  {ct+y  [Qs0+Qslv  (e)  ) sinx 


+ P Ipco+pciv  ^c^COST  + °^P 


where  p(e)  is  as  given  in  (32)  and 


p(e)v(e)  = ~[  (l-t+§-)-e'e]  (41) 

e J 

2 

The  constants  P . , P . , Q . and  Q . depend  only  on  w,  a and  w 
sj  cj  sj  c3 

and  will  not  be  listed  here.  Therefore,  the  effect  of  the  span- 
wise  correlation  is  completely  described  by  the  quantities  p(e) 
and  v(e).  Note  that  v(e)+l  as  e -*  0 . 

Having  the  steady  state  solution  for  P(x)  and  Q(x),  we 
can  now  use  (27)  to  determine  the  steady  state  meansquare  proper- 
ties of  the  blade  response.  In  view  of  (40) , a steady  state  solu- 
tion of  (27)  may  be  taken  in  the  form 


(U,S,V)  = {U,S,V}  + y [ (U  , S ,V  Jsinx 

s s s 


+ {U,S,V} cosx ] + 0(y*) 

O V C 
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The  0(1)  terns,  U,  £ and  V,  are  just  the  steady  state  solutions 

for  the.  hover  case  given  by  (33).  By  the  method  of  undetermined 

coefficients,  the  constants  U , S , ...,  V are  solutions  of  a 

s s c 

system  of  six  coupled  linear  algebraic  equations  which  nay  be 
written  as  two  sets  of  three  complex  equations 

< (1  + Ji)V  + u2U  = 2iQ  - jiV 

- (43) 

- v + (o  - j - = p + £iu 

0 - 2iS  = 0 (44) 

with 


{U , S , . . . , Q } = {US,...,Q  } - i(U  ,...fQ  } 

.2.2 


(45) 


{Ps'Pc'Qs,Qc}  “ 36AP ^Ps0' * * * 'Qc0*  + {Psl,*-*'Qcl}v(e)]  (46) 


It  follows  from  (33)  and  (45)  that  the  solution  of  (43)  and 
(44)  can  be  put  in  the  form 


Us  = a2p  (e) [Us0+Uslv  (e) ] , etc.  (47) 

2 

where  UsQ,  Uc-p'**’'vci  depend  only  on  y,  a and  w . 

In  the  case  c = 0,  the  solutions  for  the  functions  U(x), 

S(t)  and  V(t)  as  given  by  (42)  are  exactly  the  approximate 
steady  state  variances  and  covariance  of  the  flapping  response 
obtained  in  [8]  for  low  advance  ratio  flight  and  will  be  considered 
known.  Our  concern  here  is  with  the  effect  of  a finite  spanwise 
correlation  length  (c>0)  on  these  response  statistics.  This  effect 
is  completely  described  by  the  two  quantities 


p (e)  and  v (e)  . 


>» , 
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From  the  plot  of  v(e)  in  Figure  7,  we  see  that  this  monotone  de- 
creasing function  changes  by  less  than  10%  of  its  value  at  e = 0 
as  the  correlation  length  shortens  (from  infinity)  to  a fraction  of 
the  blade  length.  Therefore,  the  main  effect  of  a spanwise  load 
correlation  is  in  the  amplitude  factor  p(e).  As  both  p and  v 
decrease  with  increasing  e,  a correlation  length  shortening  in  the 
low  advance  ratio  range  gives  rise  to  a reduction  in  the  time  aver- 
age of  the  meansquare  flapping  properties  as  well  as  in  the  fluctu- 
ation about  these  average  values. 

2 

Solutions  for  O (y  ) -terms  in  (38)  and  (42)  have  also  been  ob- 
tained. In  the  e = 0 case,  these  terms  involve  the  second  harmonics 

cos  2t  and  sin  2x.  The  effect  of  a finite  spatial  correlation 

2 

length  on  these  0(y  ) -terms  is  qualitatively  similar  to  that  on  the 

2 

0(1)  and  0 ( y ) terms.  As  such  the  explicit  solutions  for  the  0(y  )- 
terms  will  not  be  given  here. 


6 . Numerical  Solution  for  Arbitrary  Advance  Ratio 


If  y is  not  small  compared  to  unity,  the  situation  is  much 
more  complicated  since  the  effect  of  reverse  flow  is  no  longer 
negligible.  For  y £ 1,  we  have  from  (20) 


c (x ) = 


|cn(T) 
lcn(t)  + 


(2m7T<T  < ( 2m+l)  7T  ) 


(48) 


(3-4cos2t+cos4t)  (( 2m+l)  r£T£  ( 2m+  2 ) tt  ) 


and 


f knh> 

5{T)  "(  4 

l k (t)  - { 2sin2i-sin4i 

n 4 o 


(2m-n<T<  (2m+l)  n ) 


) ( ( 2m+l ) it  <t  < ( 2m+2 ) tt  ) 


(49) 


Evidently,  the  effect  of  reverse  flow  is  negligible  in  c and  k 


if  y 3 <<  1 • 


From  (25)  , we  get 


r (x,t) 


l n 

t r 


n 


(x,  x) 
(x,  x ) 
(x,  x) 


- ^(x.x) 
+ (x, t) 


(2mtt£i£(2m+l)  tt  , CKx£l) 

( ( 2m+ 1 ) it<_T < ( 2m+  2 ) tt  , X£-y sinx  ) 
( ( 2m+  l)7i<x<(2m+2)7r,x>-ysinx) 


where 


(50) 


r£(x,x)  = ye  3{e  c (1-x) [ (2+2e+t2)+e (1+e) ysinx] 


Qe (x+ysinx ) 


(2-eysinx) } 


(51) 


ry(x,i)  = ye  3{e  e“  (2-eysinx ) -e  c (-+ysinx)  (2+eySinx)  } 

With  (50)  and  (51),  it  is  not  difficult  to  show  that  the  effect 

3 

of  reverse  flow  can  be  neglected  in  r(x,x)  if  y <<  1 at 

least  for  e <<  1 and  e >>  1. 

For  y > 1,  the  entire  blade  is  subject  to  reverse  flow  in 

the  range  -sinx  > — so  that 

y 

c ( x ) = -cn(x),  k(x)  = -kn(T)/  r(x,x)  = -rn(x,x)  (52) 

for  all  t in  the  range  y " t £ t £ y + 'I'  where  = cos-3'  (1/y ) . 

Having  the  expressions  for  c,  k and  r,  we  can  now  solve  the 

initial  value  problem,  (23)  and  (24) . numerically  using  a 4th 

order  Runge-Kutta  scheme  for  x,  say  Xg  = 0 ,x^  , . • • ,xm  = 1.  With 

f . (x.  ) = f ( x,  , x • ) , the  set  of  solutions  {p  . (x,  ) ,q  . (x,  ) } for  a 
jk  xj  3 X J K 

fixed  j is  used  in  (21)  to  get  P(x^)  and  Q(tj)  with  the 
integrals  evaluated  bv  Simpson's  rule.  Once  P(x^)  and  Q(x^) 
are  calculated,  the  initial  value  problem  (26)  and  (27)  is  solved 
numerically  again  by  a 4th  order  Runge-Kutta  scheme.  Within  the 
stability  boundaries  of  the  two.  sets  of  equations,  (23)  and  (27) , 


-19- 


I 


we  yet  accurate  steady  state  periodic  solutions  of  the  i.ieansquare 

blade  flapping  properties  after  four  blade  revolutions  for  the 

realistic  range  of  values  of  y (2£y£l2) . For  a fixed  set  of 

2 

y,  p , e,  a and  oj  , the  entire  solution  process  for  P,  Q,  U,  S 
and  V consumes  about  50  seconds  on  a UNIVAC  1106  if  21  stations 
along  the  blade  span  are  used  in  the  numerical  evaluation  of  the 

integrals  on  the  right  side  of  (21) . 

2 

With  = a (a  constant),  the  class  of  random 

functions  characterized  by  (3)  seems  to  adequately  describe  the 
random  inflow  associated  with  atmospheric  turbulence  at  altitude 
higher  than  300  ft.  above  terrain  if  the  effect  of  the  spatial 
variation  of  the  vertical  turbulence  component,  of  the  longitu- 
dinal turbulence  component  itself  and  of  the  blade  motion  are 
all  neglected  (see  [4]  and  references  therein).  In  that  case, 
we  have  a = 2p£,/L  where  i is  the  blade  length  and  L/2  is  the 
scale  length  of  the  vertical  turbulence  component.  L is  about 
400  ft.  for  an  altitude  of  300-700  ft.  above  terrain  and  is 
several  thousand  feet  for  higher  altitudes.  From  the  expression 
for  a,  we  see  that,  at  the  low  advance  ratio  range,  the  correla- 
tion time  is  long  compared  to  one  blade  revolution  for  existing 
blades  which  range  from  33  ft.  to  100  ft.  As  such,  the  results 
of  section  (5)  for  the  low  advance  ratio  range  serve  only  to 
indicate  the  qualitative  effect  of  a spatially  correlated  inflow; 
we  are  mainly  interested  in  the  case  of  high  advance  ratio  flight. 

The  meansquare  flapping  response  of  the  blade  to  a zero  mean 
X(x,t)  with  a correlation  function  given  by  (3)  have  been  studied 
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with  the  help  of  the  numerical  solution  scheme  outlined  in  this 

section  for  a wide  range  of  the  blade  and  load  paramenters. 

The  numerical  solution  shows  that  the  perturbation  solution  of 

2 

section  (5)  (including  0(y  ) -terms)  gives  a very  good  approxima- 
tion of  the  exact  solution  for  y ^ 0.4.  It  also  shows  that  the 
effect  of  a finite  e for  all  0 £ y <_  1.6  is  qualitatively 

similar  to  that  indicated  by  the  perturbation  solution.  The 

2 • 2 

actual  distributions  of  the  steady  state  <$  > and  <<J>  > are 

given  in  Figures  (3),  (4),  (5)  and  (6)  for  y = 1.6  and  y = 1.2 

and  for  the  two  extreme  rotor  disc  sizes,  £ = 33^-  ft.  and 

£ = 100  ft,  operating  at  300  ft.  - 700  ft.  above  terrain 

2 

(L  = 400  ft.).  We  have  taken  w = 1 in  these  examples  since 

most  existing  blades  are  hinged  at  the  blade  root.  We  see  that, 

aside  from  an  increase  in  the  magnitude  of  the  meansquare  response, 

2 *2 

an  increase  in  y tends  to  shift  the  time  when  <(j>  > and  <$  > 

attain  their  maximum  values  further  toward  the  midway  point  and 

the  end  of  the  backstroke,  respectively. 

Finally,  we  show  in  Tables  1 and  2 the  effect  of  the  Lock 

2 *2 

number  y on  the  peak  values  of  <4>  > and  <4>  >.  We  see  in 
particular  that  the  amplitude  growth  with  y is  nonlinear  and  the 
growth  rate  depends  significantly  on  y but  not  at  all  on  e. 

7.  Autocorrelation  Functions 


Having  determined  U,  S and  V,  we  can  now  calculate  the 
autocorrelation  of  the  flapping  angle  <P(t)  which  characterizes 
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the  second  order  statistics  of  the  flapping  response.  We  begin 
by  multiplying  (1)  by  w{x,/t')  and  ensemble-averaging  the  result 
to  get 

Rtt  + Y0  |x+dsinT  |rt  + LxxtRl  = Y0  I x+ysim  | A (x,x  ;x'  ,t  ' ) (53) 

where  R(x,t;x',t')  = <w (x, r ) w (x ' ,x ' ) > and 

A (x, x ; x ' , t 1 ) = <A (x, x ) w (x' , x ' ) > . To  get  the  yet  unknown  load- 
response  correlation  A,  we  multiply  (4)  through  by  w(x',x') 
and  ensemble  average  giving  us 

A + aA  = /2a  <n (x, x ) w (x ' , x 1 ) > =0  (t  > x')  (54) 

x 

where  the  right  hand  side  vanishes  for  x > x'  by  (5).  At 
x = x',  we  have  from  the  relevant  definitions 

A (x,x ' ;x' ,x ' ) = p(x,x',x')  (55) 

It  follows  from  (54),  (55)  and  the  assumption  of  rigid  flapping 

that 


(56) 


(57) 


A (x, x ; x ' ) = p (x, x ' ) e 


(x  > x') 


(58) 


-22- 


Upon  introducing  the  rigid  flapping  assumption  into  (53) , 
we  get 


Rtt  + c(t)Rt  + [w  +k(T))R  = A (t  , t 1 ) (t  > t') 


(59) 


where  R(t;t')  = <$(t)$(t')>  and 

A(t;t')  = 3yn  [ x | x+usinx | A (x, t ; t ' ) dx 
u J0 


= 3yQe  a^T  T x|x+ysinx |p(x,T 1 )dx 

' 0 


= 1 -a(T-T')- 


2 Ye 


P ( T ; T ' ) 


(60) 


with  P(x';x')  = P(t').  Note  that 

R(T';T')  = <<J>2  (T  1 ) > = U(T'),  RT(T';T')  = S(T')  (61) 

Rtt,  ( t ' ; t ' ) = V(t')  (62) 

The  two  conditions  in  (61)  serve  as  initial  conditions  for  (59) . 
But  even  without  solving  the  initial  value  problem  (59)  and  (61) 
explicitly  for  R(t;t'),  the  following  informative  observation 
can  be  made.  Since  the  effect  of  a spatial  load  correlation 
appears  only  in  P(t;t')  which  is  a periodic  function  of  t 
and  x'  at  steady  state,  the  correlation  time  of  the  response 
depends  only  on  the  parameters  a and  y and  not  on  e.  Within 
the  framework  of  rigid  flapping,  our  particular  type  of  spanwise 
load  correlation  only  modifies  the  amplitude  of  the  autocorrela- 
tion of  the  response. 


L 
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